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Abstract. Counting periodic orbits of endomorphisms on the 2-torus is considered, with 
special focus on the relation between global and local aspects and between the dynamical 
zeta function on the torus and its analogue on finite lattices. The situation on the lattices, 
up to local conjugacy, is completely determined by the determinant, the trace and a third 
invariant of the matrix defining the toral endomorphism. 



1. Introduction 

The iteration of a continuous mapping T of a compact space J7 into itself provides an 
example of a Z-action on Q and an important discrete dynamical system, usually written as 
(i7, T). When is a metric space, the system {f2, T) is called chaotic when the periodic orbits 
of T are dense in f2 and when also a dense orbit exists, see p] for details. In general, significant 
information about T is contained in the periodic orbits of T and in their distribution over 
i?. Knowledge of the periodic orbits can be used to detect characteristic properties of T. For 
example, if T' represents another continuous mapping of f2, then a necessary condition for T 
and T' to be topologically conjugate is that they share the same number of periodic points 
of each period (presuming these numbers are finite). 

It is the aim of this paper to contribute to the structure of periodic orbits and related 
issues of conjugacy for the case of endomorphisms of the 2-torus, represented as usual by the 
action (mod 1) of an integer matrix M G Mat(2,Z) on ~ M^/Z^. A well-studied subclass 
consists of the toral automorphisms, represented by elements of the group GL(2, Z), being the 
subgroup of matrices with determinant ±1 within the ring Mat(2, Z). Particularly important 
are the hyperbolic ones (meaning that no eigenvalue is on the unit circle), which are often 
called cat maps. Since these are expansive, all periodic point counts are finite [4H Thm. 5.26]. 
Hyperbolic toral automorphisms are also topologically mixing and intrinsically ergodic, see 
|23i I41j . By the Bowen-Sinai theorem, compare [15, Thm. 2.2], this has the consequence that 
the integral of a continuous function over T'^ equals its average value over the points fixed by 
M"^ in the limit as m ^ oo. 

The topological entropy of a hyperbolic toral automorphism M £ GL(2, Z) is given by 
loglAmaxI, where Amax is the eigenvalue of M with modulus > 1. This is also the metric 
(or Kolmogorov-Sinai) entropy of M, and completely determines the dynamics up to metric 
isomorphism, compare [2j. This does not imply topological conjugacy though, and one im- 
portant difference emerges from the periodic orbits, which live on a set of measure 0. Indeed, 
on T'^, it is well-known that the periodic orbits of hyperbolic linear endomorphisms lie on the 
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invariant lattices given by the sets of rational points with a given denominator n, also known 
as n-division points (see Section [2.21 for more). One of our main themes in this paper is the 
interplay between the periodic orbit statistics on a certain lattice (which we call local statis- 
tics) versus periodic orbit statistics on the union of all lattices (which we call global statistics). 
What determines when two cat maps have the same global statistics? What determines when 
two cat maps have the same local statistics on a certain lattice or on all lattices? 

At the outset, it is worth saying that there have been many investigations by others into 
classifying the periodic orbits of cat maps, spread over a diverse range of the mathematics and 
physics literature, compare [m [321 [Ml nil m linfl and further references given there. One 
motivation for this has come from the interest in spatial discretisations of dynamical systems, 
itself motivated by computer (screen) realisations of continuous phase spaces. The time of 
recurrence of a hyperbolic M £ GL(2, Z) on the toral rational lattice with denominator n is 
denoted by per(M, n), where this is the least common multiple of the periods present on the 
n-division points. The dependence of per(M, n) on n and related lower and upper bounds 
have been addressed in [T71 [THl [SH [26l [I8l [37]. In particular, the surprisingly low value of 
per(Af, n) for some high values of n (which correspond to a very fine rational discretisation 
of the torus) has been investigated in |17[ [TO] where it is shown that per (M, n) < 3n (see j37) 
for refinements of this bound) . 

A strong motivation for studying cats maps on lattices comes from quantum mechanics, 
compare [211 IMl [25l [271 [2H1 HSl [181 [16]. As described in these references, quantum cat maps 
and their perturbations are built from (classical) cat maps and their perturbations restricted 
to a particular rational lattice (called the Wigner lattice in this instance). For this reason, 
properties of a cat map that manifest themselves only on some rational lattices, but not on 
others, can induce properties of the corresponding quantum cat map on some Wigner lattices, 
but again not on others. Important cases of this occur for symmetries or (time) reversing 
symmetries of a cat map, these being automorphisms of the torus that commute with the cat 
map, respectively conjugate it into its inverse. By way of illustration, it was shown in |7] that 
the first hyperbolic toral automorphism A E SL(2, Z) which is not conjugate to its inverse 
in GL(2,Z) (which actually also excludes topological conjugacy to its inverse, compare [I]) 
occurs for trace 20. However, the global absence of time-reversal symmetry did not affect the 
statistics of the eigenvalues of the quantum cat map built from this example [25] . As explained 
there, this phenomenon is due to the fact that the quantum cat map retains (time) reversing 
symmetry because A is conjugate to its inverse mod n for any n and that the quantum cat map 
is constructed from the reduction of A mod n. Significantly, the conjugating matrix on each 
lattice depends on n, consistent with there being no global reversor. Recently, there has been 
quite some interest in dealing with this challenge of so-called pseudo-symmetries of quantum 
maps that are not quantisations of symmetries of the cat map on the torus, but instead are 
manifestations of local symmetries of the cat map restricted to some lattice \25\ [271 116] . 



These investigations have used a variety of techniques. One is tempted to say, corrupting a proverb used 
by Mark Twain and others (http://www.worldwidewords.org/qa/qa-morl.htm): There's more than one way 
to skin a cat (map). 
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When trying to sort cat maps by their global or local periodic orbit statistics, conjugacy is 
highly relevant. Conjugacy of GL(2,Z) matrices is another topic that has arisen in a broad 
variety of contexts and has been considered by many (see \40 \ 134 1 13] and references therein). 
Conjugacy is determined by a triple of invariants, namely the determinant, the trace and 
one other invariant which can be related to ideal classes, representation by binary quadratic 
forms or topological properties p5j. Conjugacy in GL(2,Z) can also be completely decided by 
using the amalgamated free product structure of PGL(2, Z), which attaches a finite sequence 
of integers to each element which corresponds to its normal form as a word in the generators 
of the amalgamated free product ^ . Clearly, COnjU^QitG Celt maps share both the same global 
period statistics and the same local statistics on each rational lattice (where the dynamics 
is conjugate via the localisation mod n of the global conjugating matrix). Also, cat maps 
that are just conjugate on a given rational lattice will share the same local statistics on that 
lattice. Being able to decide global and local conjugacies is thus clearly important, as the 
statistics is the same for all elements of a conjugacy class. But if two cat maps share the 
same local statistics on a given rational lattice, are they linearly conjugate on that lattice and 
what determines this? 

The results of this paper will go some way towards answering the questions raised above. 
After recalling some well-known facts in Section [21 we look at periodic orbit counts for integer 
matrices in terms of zeta functions in Section [3l The dynamical zeta function for the global 
counts is described by Proposition [H generalising a result of p5] . We then discuss a zeta 
function for the local periodic counts derived from the action of an integer matrix on a 
rational lattice. Theorem [1] relates the global and local zeta functions in a suitable limit. This 
is followed by an interpretation in terms of finite Abelian groups. Section [4] addresses the 
issue of local conjugacies of linear endomorphisms on rational lattices. The matrix gcd, which 
we define in Section BTl turns out to be a key quantity. It is preserved by GL(2, Z) conjugacy, 
so it provides a quick tool to see that two GL(2, Z) matrices with different matrix gcd are not 
conjugate on the torus. On the other hand, Theorem[2]and Corollary [3] show that two integer 
matrices that share the same determinant, trace and matrix gcd are linearly conjugate on all 
rational lattices of the torus. As an illustration of this result, consider our discussion above 
of quantum cat maps and time-reversal symmetry. The fact that any M £ SL(2, Z) shares 
determinant, trace and matrix gcd with means that the two matrices are conjugate on 
all rational lattices, though not necessarily by matrices that derive from one and the same 
matrix on the torus. This is nevertheless sufficient to guarantee that the associated quantum 
cat map has time reversal symmetry. 

2. Notation and general setting 

Here, we describe our setting and recall some well-known facts, tailored to the situation at 
hand. While we go along, we also introduce our notation and establish further connections 
with related topics in the recent literature. 

2.1. Counting orbits. Consider a compact space f2 and some (continuous) mapping T of 
into itself. Let Fixm(r) := {x € X \ T^x = x} be the set of fixed points of T™. Of particular 
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interest are the fixed point counts, defined as 



(1) 



card{2; G n \ T"^x = x} = card(Fix„(r)) 



which need not be finite in general. In many interesting cases, including all expansive home- 
omorphisms, this is the case though, including the toral endomorphisms without eigenvalues 
on the unit circle. 

The quantity am has the disadvantage that one keeps recounting the contributions for 
all i\m. Clearly, the fixed points of genuine order m permit a partition into disjoint cycles, 
each of length m. If Cm is the number of such cycles, one thus has the relation 



d\Tn 

An application of a standard inclusion-exclusion argument, here by means of the Mobius 
inversion formula from elementary number theory, results in the converse identity. 



where fj,{k) is the Mobius function, compare \22\ [33] and references therein for details. 

Remark 1. Recall from [33j that a sequence (am)mGN of non-negative integers is termed 
exactly realised when it is the sequence of fixed point counts of a (continuous) dynamical 
system. This happens if and only if the derived sequence {cm)meN is a sequence of non-negative 
integers [331 Lemma 2.1]; see [42j for interesting examples other than toral endomorphisms 
and |3iJ for recent extensions of the concept. 

For later use, we briefly summarise some properties of the fixed point and orbit counts. 
Let (am)meN and (cm)mGN be a matching pair of such sequences, hence related as in Eqs. ([2]) 
and ([3D. The sequence of fixed point counts is called periodic, when an n G N exists so that 
am+n = cLm holds for all m G N. The least n with this property is called the period of the 
sequence (am)meN- The following consequences are standard. 

Fact 1. Let the non-negative integer sequences (am)mgN cind {cm)meN satisfy Eq. ([2]). If am 
is periodic with period n G N, one has Cm = for all m > n. Conversely, if only finitely many 
orbit counts Cm differ from 0, am is periodic, with period n = lcm{m G N | 7^ 0}. □ 

In extension of the usual practice for automorphisms, compare [23], we call a toral endo- 
morphism M G Mat(2,Z) hyperbolic when it has no eigenvalue on the unit circle Recall 
that the standard 2-torus is ~ R^/Z^, where is the square lattice in the plane. It is a 
compact Abelian group, which can be written as := [0, 1)'^, with addition defined mod 1. 

Fact 2. IfMe Mat(2,Z) is hyperbolic, FiXm(M) C (T^n Q^) holds for all m G N. 

Proof. This can be shown by the concrete argument used in [23[ Sec. 1.8, p. 44]. In essence, 
for any m G N, the equation M^x = x mod 1 means (M"^ — l)x = v for some integer vector 
V. Since 1 is not an eigenvalue of M™, the integer matrix M™ — 1 is invertible, with rational 
inverse. Solving for x then gives the claim. □ 



(2) 




(3) 




dim 
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This property motivates to look at periodic points of toral endomorphisms on certain 
subsets of n as well. 

2.2. Invariant lattices on the 2-torus. Recall that a lattice in a locally compact Abelian 
group is a co-compact discrete subgroup, such as in M^. For T^, a lattice simply means a 
discrete subgroup of it, which is then a finite Abelian group. The relevant lattices on are 

(4) := {{^,^)\0<k,£<n} , forn G N, 

also known as n-division points, because they are invariant under toral endomorphisms. 

Consider an integer matrix M that acts on (meaning that it acts via matrix multipli- 
cation, evaluated mod 1). Clearly, one then has ML„ C L„, and interesting information on 
the orbit structure of the toral endomorphism or automorphism (det(M) = ±1) is contained 
in the distribution of its orbits on L„. Alternatively, one can characterise -L„ via 

(5) Ln = {xeT'^\nx = (mod 1)}. 

Remark 2. When looking at the action of M on Ln numerically, it is usually easier to replace 
Ln by Ln '■= {(^,1) I < k,i < n}, with the equivalent action of M defined mod n. This also 
applies to various theoretical arguments involving modular arithmetic. Consequently, we use 
Ln (with action of M mod 1) and L„ (with action mod n) in parallel. 

From now on, we use the abbreviation Z„ = TLjnL for the finite integer ring mod n, and 
= {1 < A; < n I gcd(/c,n) = 1} for its unit group. Let now n > 2 be arbitrary, but fixed. 
If we have a matrix A E GL(2,Z), its reduction mod n is still invertible in Mat(2,Z„), and 
thus an element of Mat(2,Z„)^, the unit group within Mat(2,Z„). This group is often also 
called GL(2,Z„), though its elements need not have determinant ±1. 

Remark 3. A matrix M G Mat(2,Z) that is not a toral automorphism (so det(M) ^ ±1) 
may still be invertible on a given lattice L„ (or on L„), meaning that M mod n is an element 
of Mat(2,Z„)^. This happens precisely when det(M) G Z is relatively prime to n, which is 
equivalent to det(M) G Z^ (where det(M) is taken mod n). 

When M is a toral automorphism, or a toral endomorphism with det(M) G Z^, the 
invertible action of M on the finite set L„ induces a "permutation of the r? elements of -L„, 
with one element (x = 0) always fixed. So, the induced permutation vr^^ can either be viewed 
as an element of the symmetric group 5^2 or of 5'„2_i. The permutation is of finite order, 
which must divide (n^ — 1)' by Lagrange's theorem. The actual order of M on L„ is given by 

(6) ord(M,n) := gcd{m G Nq | = 1 mod n\. 

Clearly, ord(M, 1) = 1 in this setting. When M is not invertible on L„, the definition results 
in ord(M, n) = 0; otherwise, ord(M, n) is the smallest m G N with = 1 mod n. 

By an application of Fact [H the orbit count sequence on L„ has vanishing entries beyond 
index m = ord(M, n). It is natural to define af^ := cardja; G L„ | M^x = x (mod 1)} and 

c^"^ := cardjcycles of M on L„ of length £} 
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as the induced fixed point and orbit counts on again mutually related as in Eqs. ([2]) and 
([3]). It is not hard to see that 

(7) ord(M,n) = per(M,n) := lcm{m g N | c^^ 7^ 0}, 

meaning that per(M, n) is the 1cm of the lengths of the cycles on L„. Clearly, we can only 
have Cm^ 7^ when m is a divisor of per(M, n). 

When M is a toral endomorphism with gcd(det(M), n) 7^ 1, it is not invertible on the 
lattice Ln- Consequently, not all points of -L„ show up in periodic orbits of M, and one 
has 'pretails' to the periodic orbits (there is always at least one periodic orbit on L„). Two 
matrices on L„ may thus share the same cycle structure there, but show different pretail 
patterns. As the latter give rise to a directed pseudo-graph on Ln (where 'pseudo' simply 
refers to the possibility of directed loops at a vertex), with the points of Ln as vertices and the 
directed edges derived from the action of M, it is reasonable to coin an adequate definition 
that covers both the case when M is invertible on L„ and when it is not. 

Definition 1. We say that two matrices M and M' have the same local statistics on L„ when 
the induced directed pseudo-graphs on L„ are isomorphic as graphs. When M and M' are 
invertible on Ln, this is equivalent to saying that they have the same periodic orbit counts. 
Otherwise, this is equivalent to saying that they have the same periodic orbit counts plus 
isomorphic pretails between corresponding periodic orbits of the same length. 

This definition shows that two matrices with different pretail structures for equal orbit 
counts would not have the same local statistics. Note that, when M and M' are invertible on 
Ln, they have the same local statistics if and only if their associated permutations vr^"^ and 
vr^"^, are conjugate in 5„2. 

2.3. Powers of 2 x 2-matrices. Let M G Mat(2,C) be a non-singular matrix, with D := 
det(M) 7^ and T := tr(M). Define a two-sided sequence of (possibly complex) numbers Pm 
by the initial conditions p_i = —1/D and = together with the recursion 

Pm+i = Tpm - Dpm-1, for m > 0, 

(8) 1 

Pm-1 = J^{Tpm-pm+l), for m < -1. 

This way, as D 7^ 0, pm is uniquely defined for all m £ Z. Note that the sequence {pm)mez 
depends only on the determinant and the trace of M. When M £ Mat(2, Z), one has Pm S Q, 
and Pm G ^ for m > 0. When M G GL(2, Z), all Pm are integers. 

Let us first note an interesting property, which follows from a straight-forward induction 
argument (in two directions). 

Fact 3. The two-sided sequence of rational numbers defined by the recursion ([8]) satisfies the 
relation p^ — Pm_^iPni_i = D^~^ , for all m G Z. □ 

Lemma 1. Let M G Mat(2,C) be non-singular. For m G Z, one has the relation 

= PraM-Dpm^il. 

In particular, one has = j){Tl — M). 
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Proof. The initial conditions p^^ = —1/D and Pq = imply that the relation = 1 is 
matched. Let us first look at the positive powers of M. Assuming the claim to hold for some 
integer m > 0, we can use the Cayley-Hamilton theorem to proceed inductively. Indeed, from 
M2 = tr(M) M - det(M) 1 = TM -Dl, one finds with §^ that 

= {Tpm - Dpm-l)M - Dpml = Pm+lM - Dpml, 

which establishes the claim for all m > 0. The special relation for is just a reformulation 
of the inverse of a 2 x 2-matrix, while the statement about the negative powers follows from 
another induction argument, using the reverse recursion for the Pm with negative index. □ 

Remark 4. When M £ Mat(2, C) is singular, so D = det(M) = 0, one can still meaningfully 
define the numbers pm for all m > 1. In fact, they are then simply given by pm = T™~^, 
with T = tr(M). The formula from Lemma [T] simplifies to M'" = PmM = T^~^ M, which is 
valid for all m > 1, while Fact [3] remains true for all m >2. The numbers pm are particularly 
useful to determine the periods ord(M, n) of Eq. ([6]). 

3. Dynamical zeta functions and periodic orbit statistics 

To deal with combinatorial quantities such as the fixed point counts Om, it is advantageous 
to employ generating functions. They provide a nice encapsulation of these numbers and 
permit the derivation of several asymptotic properties as well. Here, the concept of a dynam- 
ical zeta function, compare [36J, is usually most appropriate. Consequently, given a matrix 
M £ Mat(2,Z), we set 

oo 

(9) CmW := exp(^^i-), 

where, from now on, := cardjx G FiXm,(M) | x is isolated} is the number of isolated fixed 
points of M"^. We say more about this below. 

Remark 5. The ordinary power series generating function for the counts Om can be calculated 
from Ca/(0 Ylm>i '^rn.t"^ = ^ ^ log (Cm (^)) • '^^^ significance of the formulation used in 
Eq. ([9]) follows from the fact that it has a unique Euler product decomposition [36j as 

(10) At) = ^ (^-*'") = n (i-n^^ 

' cycles C m>l 

where \C\ stands for the length of the cycle C and Cm is now the number of isolated cycles of 
M on of length m, as determined from Formula dSD- Consequently, the role of cycles in 
dynamics is similar to that of primes in elementary number theory. 

3.1. Global considerations. Let M G Mat(2,Z) and observe that Fixm(M) is an Abelian 
group. In fact, it is a closed subgroup of T^. Consequently, when 1 happens to be an eigenvalue 
of M"^ for some m, there is a continuous subgroup of fixed points of M™, and one cannot 
have any isolated fixed points in addition, due to the group structure of FiXm(M), see [6] for 
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details. For a given toral endmorphism, any set FiXm(M) is thus either finite (when is an 
isolated fixed point) or a continuous submanifold of T^. 

With hindsight, this motivates our definition in ([9]), and we now need to have an explicit 
formula for the number Um of isolated fixed points of M G Mat(2, Z) on T-^. This is possible via 
a standard argument that involves areas of fundamental domains of planar lattices, compare 
[H [15] . The subtorus case is treated in [6l Appendix] . 

Fact 4. If M £ Mat(2,Z) is an arbitrary integer matrix, the number of isolated fixed points 
of M"* on is given by 

am = |det(M™-l)|, 

which is valid for all m € N. In particular, all these counts am and the corresponding cycle 
numbers Cm are finite. Moreover, whenever Om = 0, no isolated fixed points exist, and one 
has subtori of fixed points instead. □ 

One CcLii express Qjfji in terms of the numbers p^ji from dH). Indeed, when D = det(M) / 0, 
by inserting the formula from Lemma [1] and by also using Fact El one finds 

(11) det(M"^ - 1) = -pm+l + Dpm-i + i?™ + 1, 

which is valid for all m > 1. By Remark HI one can check that inserting D = gives the 
correct formula also for det(M) = 0, namely det(M'^ - 1) = 1 - T™, with T = tr(M). With 
these formulae, one can determine the dynamical zeta functions of the isolated fixed points 
for many cases explicitly, compare [39, Sec. 1.4] for an alternative approach. 

Example 1. Let us start with M E Mat(2,Z) being singular, and of trace T. One finds 
(^j^^(t) = 1/(1 - t) for T = and 

(12) CM(t) = 

for r 7^ (note the special role of T = 1, with Cm(^) ~ existence of subtori 

of solutions). The explicit derivation follows from Remark |1] and the standard power series 
identity (with ^3 = 1 as radius of convergence) 



(13) log(l -z) = -Y, 



, m 

m=l 



When M = kl with k € 7j, a simple calculation results in 

which is also valid for M = (9^)1 with 6 E Z arbitrary. If A; = 1, we are back to a case with 
subtori of fixed points, again reflected by Ca/I*) = ^■ 

More generally, consider a non-singular upper triangular matrix of the form M = ( q ^ ) 
with ad 7^ 0, so that 6 = sgn(L') 7^ 0. Then, one finds the zeta function 

detjl -6tM) _ {l-6at){l-6dt) 
^ ^ " {l-5t){l-5Dt) ~ {l-6t)il~5Dt)' 

again using standard manipulations of power series. 
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3.2. Zeta functions of toral automorphisms. For M G GL(2, Z), Fact U] permits a deriva- 
tion of the dynamical zeta function as follows, a special case of which was also given in [15j . 

Proposition 1. Let M G GL(2,Z) he hyperbolic, and define a = sgn(tr(M)). Then, with 
the coefficients am = cardjx G | M™x = x (mod 1)}, the dynamical zeta function Q of 
M on is given by 

{l-at){l-at det(M)) _ (1 - at)(l - a det(M) t) 
~ det(l - (jtM) ~ l-|tr(M)|t + det(M)t2- 

In particular, Cm(0 ^ rational function, with numerator and denominator in 7L^. The 
denominator is a quadratic polynomial that is irreducible over Z. Its zero tmin closest to 
gives the radius of convergence of (^^.^(t), as a power series around 0, via r^ = |tn 

Proof. Recall first from [36j that, for arbitrary matrices A G Mat(n,C), 



Tnm 



oo 

(16) exp(^ 



k J det(l - tA) ' 

with convergence for \t\ < l/g, where g is the spectral radius of A. 

If M is hyperbolic, the general formula for the am from Fact U] can be used to derive 

am = cj'"(tr(M'") - (1 + det(M)™)) 

by observing that the two eigenvalues of A can be written as A and det{A)/\. For the detailed 
argument, one may assume |A| > 1 and check the different cases. Note that a hyperbolic toral 
automorphism is never of trace 0. 

The formula for the zeta function now follows from by inserting the expression for am 
and using the relation (jl6p together with the power series identity (jlSp . The statement on 
the nature of the rational function is then clear. With M G GL(2,Z), the denominator only 
factorises for tr(M) = 0, det(M) = -1 or for tr(M) = ±2, det(M) = 1, both cases being 
impossible for hyperbolic matrices. The result on the radius of convergence is standard. □ 

Example 2. Probably the best known hyperbolic toral automorphism is the one induced by 
the 'classic' or golden cat map 

^-(? 

It has det(M) = —1 and is thus orientation reversing (sometimes, as in [23j, its square is used 
instead). Prom Proposition [1] or from [6], one obtains 

l-t2 



m>l 



with am = fm+1 + fm-1 - (l + (-1)") and Cm according to Eq. see also entries A001350 
and A060280 of [35]. Here, fm are the Fibonacci numbers, defined by the recursion fm+i = 
fm + /m-i) for m > 0, together with the initial condition /o = and /_i = 1. The first few 
terms of the counts are given in Tabled) As an aside, note that Ca/(0 = 1 + Ylm=of^^"^^ 
and one has M™ = fm M + fm-i !> the latter being valid for all m G Z. 
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m 


1 


2 


3 


4 


5 


6 


7 


8 9 


10 


11 


12 


13 


14 


15 




1 


1 


4 


5 


11 


16 


29 


45 76 


121 


199 


320 


521 


841 


1364 


C-m 


1 





1 


1 


2 


2 


4 


5 8 


11 


18 


25 


40 


58 


90 



Table 1. Fixed point and orbit counts for the golden cat map. 

Remark 6. In Example [2l when taken mod n, the powers M"^ are periodic in m with period 
Pn, which are known as the Pisano periods, see entry A001175 of [38j. The periods of the 
sequences (am'')^g|^ divide P„, but can be smaller, as happens for n = 4 (with 3 versus 
P4 = 6) or for n = 5 (with 4 versus P5 = 20). 

Let us return to our general discussion. From Proposition [H it is clear that two hyperbolic 
GL(2, Z)-matrices with the same trace and determinant possess the same dynamical zeta 
function, hence the same fixed point counts. The converse is slightly more subtle. 

Corollary 1. Let M,M' £ GL(2,Z) represent two hyperbolic toral automorphisms which 
have the same fixed point counts. Then, Cm(^) ~ Cm'(^)- This implies det(M') = det(M) and 
either tr(M') = tr(Af) or tr(M') = -tr(M), the latter together with det(M) = -1. 

Proof. The claimed equality of the zeta functions is clear. As the denominator is irreducible 
over Z by Proposition [H we get det(M') = det(M) and |tr(M')| = |tr(M)|. Equating the 
numerators results in the two possibilities stated. □ 

Remark 7. The second possibility of Corollary [T] is realised by any orientation reversing 
GL(2, Z)-matrix M together with M' = —M. One can also check this property by an explicit 
calculation, using Fact [4] in conjunction with Lemma [Hand Fact [3l 

A concrete example is provided by the golden cat map (or Fibonacci matrix) M = ( 5 } ) of 
Example [2j Since —1 is a lattice symmetry for all L„, it is clear that also the local statistics 
of M and —M is the same on all lattices. 

The formula of Proposition [T] does not necessarily hold for the elliptic and parabolic ele- 
ments of GL(2, Z). However, those cases can be derived from the formulas given in ExamplelU 
and result, in our setting, in the dynamical zeta function for the isolated fixed points. 

3.3. Generating functions on lattices. Let us now consider a toral endomorphism M on 
the lattice L„, for some n € N. It is clear that is mapped onto itself under the action 
of M (mod 1). Alternatively, by Remark [2l we may consider Z/„, and thus the action of the 
reduction of M mod n. Recalling that card(-L„) = card(Z/„) = n^, the following observation 
is obvious. 

Fact 5. Let M £ Mat(2,Z) and n £N be fixed. Then, only finitely many of the orbit counts 

(n) 

on Ln can be non-zero, and the dynamical zeta function for the action of M on Ln has 
the property that 



is a finite product and defines a polynomial in Z[t] of degree at most r? . □ 
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Note that the degree of Zn{t) can indeed be less than nP. This happens whenever the 
action of M is non-invertible on which manifests itself in the existence of pretails to the 
actual orbits. Otherwise, Fact [5] can be sharpened as follows. 

Proposition 2. Let M £ GL(2,Z) or, more generally, let det(M) for M G Mat(2,Z) be 
coprime with n, meaning that the reduction of M is an element of Mat(2,Z„)^ = GL(2,Z„). 
Then, the dynamical zeta function (t) satisfies 

1 = n {i-*'r'"'=:z„(*). 

£ I per (M,n) 

where Zn{t) E Zi[t] has degree . In particular, with N = per(M, n), one has 



. — a AT — ri 



and the minimal period of the sequence ('^m^)^gj^ is a divisor of N . □ 

Concerning the last statement, it is worthwhile to mention that the minimal period of the 
fixed point counts on L„ can actually be smaller than per(M, n), as we saw in Remark [6j 
Note that (1 — t)\Zn{t) for all n > 1, since x = is a fixed point of M on all L„. Moreover, 
771 1 n implies Z^a \ Zn . This induces a partial order on the polynomials Zn, which permits us to 
consider the direct (or inductive) limit of them, considered within the ring of formal power 
series, compare [201 Sec. 1]. In fact, one simply has 

limZ„(t) = lim lcm{Zi(t), Z2(t), . . . , 

' n— >oo 

where 1cm stands for the least common multiple, which is well-defined here because the 
polynomial ring Z[t] is factorial, see |29j . It is clear by construction that this limit must 
divide l/Cjvf (t), the latter written out as an infinite Euler product and thus as an element of 
Z[[t]]. In our setting, the formal power series have positive radius of convergence, so that one 
can also extract asymptotic properties of their coefficients by standard tools from complex 
analysis, compare [HI |20] for examples. 

Theorem 1. If M £ Mat(2,Z) is hyperbolic, one has limZ„(t) = l/Cj\/(0 • 

Proof. When M has no unimodular eigenvalue. Fact H] gives a formula for the number of all 
fixed points of Af", which is finite. It is clear that the finitely many fixed points of Af" are 
isolated. Viewing each polynomial Zn{t) as an element of Z[[t]], it is clear that 

limZ„(t) I 

where 1/C,^.^{t) is written as the Euler product of (jlOp . hence as an element of Z[[t]], which 
contains each Zn{t) as a factor. It remains to show that each factor (1 — t"^y^ of 1/Cjvf(i) 
divides for some n G N (and then also for all multiples n' of n). 

Note that all fixed points of M*" he in Q^nT^ ~ Q^/Z^ by Fact[2l Consequently, there must 
be some n = n{m) such that all these fixed points are in L„, which implies (1 — t"^)'^\Zn{t). 
The claim now follows from the general structure of the direct limit. □ 
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Remark 8. The direct limit does not exist for all matrices in Mat(2,Z). This relates to the 
observation that endomorphisms with subtori of fixed points of some order have no isolated 
fixed points of the same order, while the intersections with the lattices L„ are still finite sets. 
In this situation, Cm(^) encapsulates the isolated fixed point counts only (if any), while Zn{t) 
gradually explores also the non-isolated ones. 

3.4. Group theoretic interpretation. Both the torus and its lattices L„ are compact 
Abelian groups (the latter even being finite). It is thus natural to also expect some group 
theoretic interpretation of counting orbits of a hyperbolic toral endomorphism M on these 
groups. Recall that 

(17) Am ■■= Fix„(M) = {x G I M'^x = x (mod 1)} 

is a subgroup of T^, hence a finite Abelian group of order a^n due to the assumption of 
hyperbolicity. Knowing am, however, generally tells us rather little about Am as a group. To 
improve on this, we employ the elementary divisor theorem for finite Abelian groups, compare 
O Thm. 1.4.2] or [29, Ch. 1.4.8]. To this end, consider 

A^m^ := Am n = {x £ Am \ nx = (mod 1)}, 

which defines a family of Abelian groups, with am^ = card (A^^) < Am^ C Am- In 

fact, since we assume here that M has no eigenvalues on S^, one has Am = UneN^^^- 

Proposition 3. The structure of the finite Abelian group Am of (I17p is completely determined 

(n) 

by the numbers Om with re S N. 

Proof. Fix rre G N, and write G = Am- Choose an isomorphism 

s 

i=l 

which exists by the elementary divisor theorem, with s G N; note that the primes pi need not 
be distinct. Set G{n) = Am^ and g{n) = am\ where g{l) = 1. In view of the elementary 
divisor theorem, it now suffices to show that, for each prime power p^ with r > 1, the power 
of p in 

gip^-^)g{p^+^) [Gip^+^) : Gipn] 
equals the number of indices i so that Pi = p and ii = r. This follows from g{p^) = where 

t = min(r, ii) . 

i:pi=p 

This uniquely specifies all elementary divisors. □ 

4. Global versus local conjugacy and orbit statistics 

The determinant and the trace of an integer matrix are not changed under GL(2, Q) conju- 
gacy. But for matrices M € Mat (2, Z), it has been known for a long time that the determinant 
and the trace are neither a sufficient nor a maximal set of invariants (this goes back to con- 
tributions by Latimer, MacDuffee, Taussky and Rademacher - see [121 El El] and references 
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therein). There are various ways of deciding GL(2, Z)-conjugacy, amounting to exploiting a 
third and final conjugacy invariant [H El [71 [3l]. It is clear that there are many interesting 
connections to class groups and class numbers of quadratic number fields, see [ID] for details. 

Example 3. Consider the two GL(2, Z)-matrices 

which share D = —1 and T = 9. One can check explicitly that the integral matrices X which 
satisfy MX = XM' are integral linear combinations of A and i?, where 

For each integer n > 2, one can find an X with det(X) coprime to n, so that the reductions of 
M and M' mod n are Mat(2, Z„) ^-conjugate. However, when taken mod 5, det(X) is always 
congruent to 0, 2 or 3. Consequently, no X exists with det(X) = ±1, whence M and M' are 
not conjugate in GL(2, Z). 

Obviously, two GL(2, Z)-conjugate hyperbolic toral automorphisms do possess the same 
dynamical zeta function on T^, equivalently have the same sequence of fixed point counts. 
They also have the same local statistics on all lattices L„. The latter statement follows from 
the observation that the conjugating GL(2, Z) -matrix element leaves all L„ invariant and 
induces a local conjugacy on all of them. This is a particular case of a result for endomorphisms 
(recall Remark [3] and its preceding paragraph): 

Fact 6. Let n > 2 he an integer. If two matrices M,M' £ Mat(2,Z) are Gh{2, 'Z)-conjugate, 
their reductions mod n are Mat{2, Zn)^ -conjugate. 

Proof. By assumption, we have M' = AMA~^ for some A G GL(2, Z), which mod n is turned 
into an equation of the same type within Mat(2, Z„), with A G Mat(2, Z„)^ . □ 

Conversely, if two hyperbolic matrices M, M' G Mat(2, Z) share the same fixed point counts 
on all lattices L„, they must also have the same fixed point counts on T^. If M, M' G GL(2, Z), 
Corollary[T]implies that they have the same determinant. They also have the same trace if they 
are orientation-preserving (their traces may differ in sign if they are orientation-reversing). 
Even in the orientation-preserving case, the equivalence of local statistics for all n does not 
imply GL(2, Z)-conjugacy. For instance, M and M' = must have the same set of fixed 
point counts on all lattices L„, as simply runs backwards through the orbits of M. But 
a hyperbolic toral automorphism need not be conjugate to its inverse: 

Example 4. The two matrices 

with D = 1, T = 20 are not conjugate within GL(2,Z), as one can check by an explicit 
calculation [7[ (in fact, this means they are not even topologically conjugate, see [HJ Fact 1]). 
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Note that the companion matrix C for both M and M ^ is conjugate to its inverse: 

This means that C and its inverse are in yet another GL(2,Z) conjugacy class than M and 
(one can calculate that there are altogether 5 conjugacy classes for D = 1 and T = 20, 
compare [5l Ex. 17]). We thus see that the set of integer matrices with the same local statistics 
on all lattices L„, can encompass more than one matrix conjugacy class on T^. 

Two conjugate matrices possess equivalent orbit structures, including pretails of periodic 
orbits. In view of Remark [2l the following local property is then clear. 

Fact 7. Let n £ N. When two integer matrices M and M' are Mat{2, Zn)^ -conjugate, the 
corresponding toral endomorphisms have the same local statistics on the lattice Ln, in the 
sense of DefinitionUl □ 

In the remainder of this section, we answer the question when two integer matrices 
(20) M = ']] and M' ^' 



c dj \c' d' 

from Mat(2,Z) are locally conjugate (mod n) for all n and hence possess the same local 
statistics (in the sense of Definition [T|) on all lattices L„. This will only depend on the 
determinant, the trace and a new invariant that we introduce next. 



4.1. The matrix gcd. Consider a 2x2-matrix M as in (pO|) . 

Definition 2. If M G Mat(2,Z), the quantity 

mgcd(M) := gcd(6, c, d — a), 

is called the matrix gcd of M, or mgcd for shor1@. Here, we take the gcd to be a non-negative 
integer, and set mgcd(M) = when b = c = d — a = 0. 

The last convention matches that of the ordinary gcd, and is compatible with modular 
arithmetic. The following consequence of this definition is obvious. 

Fact 8. For M £ Mat(2,Z), the following statements are equivalent: 

(a) The matrix gcd satisfies mgcd(M) = 0. 

(b) M = kl for some k e Z. 

(c) The minimal polynomial of M is of degree 1. 

Consequently, whenever mgcd(M) = r G N, M cannot be a multiple of the identity, and its 
characteristic and minimal polynomials coincide. □ 



generalisation of the mgcd to nxn integer matrices, n > 2, is the quantity m of [12], which is used to 
describe the least normal subgroup in GL(n, Z) containing a given element. 
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There are various other useful properties of the matrix gcd. It is immediate that 

(21) mgcd(-M) = mgcd(M) = mgcd(M*) 
holds for arbitrary M G Mat(2,Z). Moreover, with A; G Z, one has 

(22) mgcd(A;M) = A;mgcd(M) and mgcd(M + A;l) = mgcd(M). 

Finally, if M is invertible, its inverse is = ^^^^^^^ (_^^ ~a)- Consequently, for all matrices 
M £ GL(2,Z), one has the relation 

(23) mgcd(M"^) = mgcd(M) 



in addition to (j2ip . 

Most significantly, the matrix gcd satisfies the following invariance property, which can be 
seen as a consequence of the close relationship to the theory of binary quadratic forms [43j . 





f I \ 












'I] 




\d-a) 




\d-a 



Lemma 2. If M,M' £ Mat(2,Z) are two integer matrices that are conjugate via a GL(2,Z)- 
matrix, one has mgcd(M') = mgcd(M). In particular, the matrix gcd of Definition [E is 
constant on the conjugacy classes of GL(2,Z). 

Proof Assume M' = AMA'^ with A= ("5) £ GL(2,Z). One can check by a straight- 
forward calculation that this implies the linear equation 

c' = -7^ 5^ -Id 

^d'-a'l ' \2a7 -2(35 /?7 + 

where = N{A) is an integer matrix because det(^) = ±1. As the primed quantities are then 
integer linear combinations of the unprimed ones, gcd(6, c, d—a) divides gcd(6', c', d' — a'). It is 
easy to verify that det(A) = 1, so that N £ GL(3, Z) and is an integer matrix as well (this 
can also be seen directly from observing that M = A~^M'A and hence N{A~^) = {N{A))~^). 
Consequently, our previous argument implies that gcd(6', c', d' — a') divides gcd(6, c,d — a) as 
well. Within N, we thus obtain gcd(&', c', d' — a') = gcd(6, c,d — a) as stated. This conclusion 
also holds when one gcd (and then also the other) vanishes, adopting the usual convention 
that 0|0. The second claim is now obvious. □ 

Remark 9. Two integer matrices with different mgcd cannot be GL(2, Z)-conjugate. Note, 
however, that M, and C of Example H] all have mgcd = 1, but are in distinct GL(2,Z) 

conjugacy classes as discussed. <0 

4.2. Local conjugacies and a binary quadratic form. With a view to determining con- 
jugacies on L„, meaning conjugacy via Mat(2,Z„)^, consider the integer matrices 

(24) M = ^] and C 



c dj \l T 

with D = det(M) and T = tr(M). Here, C is the standard companion matrix for the 
characteristic polynomial 

(25) -Tx + D 
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of the matrix M. Let us assume that M is not a multiple of 1. To investigate possible 
conjugacies, let ^ = be another integer matrix. It is easy to check that 

(26) MA = AC ^ (^^ = m(^^. 
Whenever this semi-conjugacy holds, one has 

(27) det(A) = ca^ + (d - a)a7 - 67^ =: Q^iot^-i), 

which brings us in contact with the classic theory of binary quadratic forms. 

In fact, the quadratic form Qj^ is the key for a hierarchy of conjugacies. Clearly, M and 
C are GL ( 2, Q) -conjugate if Qj^^ represents some ^ 7^ 0. This is true unless mgcd(M) = 0, 
which relates back to Fact [8j Whenever Qj^ represents 1 or —1 (which then automatically 
implies a and (3 to be coprime), the matrices M and C are GL(2, Z)-conjugate (this is one 
of the approaches to GL(2, Z)-conjugacy taken in [HI?]). When Qj^^ represents some ^ E N, 
and when n G N is another integer that is relatively prime with /i, one has (j, G Z^, wherefore 
A is invertible in Mat(2,Z„), compare Remark[3l and the reductions of M and C mod n are 
Mat(2, Z„) ^-conjugate. 

The discriminant of the binary quadratic form Qj^j of (j27p is 

(28) A = {d-af + Abe = (a + df - 4{ad - be) = T^ - AD. 

The form Qj^ is called primitive when gcd(6, c, d — a) = 1, which means mgcd(M) = 1. 
Moreover, a representation k = Qj^j{a,^) is called primitive (or proper) when gcd(a,7) = 1. 
We need the following fundamental result from the theory of binary quadratic forms, compare 
[I3l Prop. 4.2]. 

Fact 9. Let n G N 6e arbitrary, but fixed. If the binary quadratic form Qjyj is primitive, it 
ean primitively represent some integer k that is relatively prime to n. □ 

Indeed, when a is the product of all primes that divide n, but not b, and 7 the product 
of all primes that divide n and 6, but not c, one sees that k = Qj^j{a,"f) is an integer that is 
relatively prime to n. 

Remark 10. When the quadratic form Qj^^ fails to be primitive, its discriminant A is the 
product of (mgcd(M))^ with the discriminant Z\' of the 'primitive part' of Q^y,/- Using (I28p . 
this relates to the following property of the eigenvalues of M, which are the roots of (j25p : 

A± = i(T± = ^(T±v^) = i(T±mgcd(M)v^) = i(T±5VZ^). 

In the last equality, A" is the square free part of A, highlighting that mgcd(M)|S'. In particu- 
lar, mgcd(M) = 1 whenZ\ itself is square free. The approach of [32j to studying the dynamics 
of GL(2, Z) matrices on Ln involves relating it to the multiplication of the eigenvalue A4. on 
an associated ideal. The latter is an order in a quadratic number field, but generally not its 
maximal order; see Ch. 2] and [13] for details on the connection between quadratic forms 
and orders in number fields. 
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4.3. Complete determination of local conjugacy. We now give a complete description 
of the local conjugacy problem, in the sense of Definition [TJ This makes use of a particular 
normal form over the ring Z„, with a sequence of elementary arguments. 

Proposition 4. Let M G Mat(2,Z) be a matrix with mgcd(M) = /i 7^ 0, and let C he the 
corresponding companion matrix, as in (I24p . Then, for all integers n > 2 that are relatively 
prime with n, the reductions of the matrices M and C mod n are 'M.&t{2,'Ln)^ -conjugate. In 
this case, M and C share the same local statistics on L„. 

Proof. With mgcd(Af ) = ^, the integer quadratic form ^Q^;^ is primitive. By Fact [U we can 
thus find a,/3 G Z with gcd(a,/5) = 1 and ^{a, (5) = k where k is an integer relatively 
prime with n. So, Qj^,j{a,l3) = kfj,, which is still relatively prime with n by assumption. 

This means that there is a matrix A £ Mat(2, Z) with det(A) = kfi whose reduction mod n 
is an element of Mat(2,Z„)^. Consequently, it defines a conjugacy of M and C on L„, again 
when taking all matrices mod n. The final claim is clear from Fact [71 □ 

Note that we need not consider the trivial lattice, Li = {0}, as the point x = is fixed by 
all matrices. 

Example 5. The conjugation in Proposition d] need not be SL(2, Z„)-conjugacy, that is, the 
conjugating element need not have determinant 1 (mod n). To illustrate this, consider the 
matrices 

with D = —2, T = 4 and, from (f27l) . quadratic form (5^(0,7) = 2a^ — 87^ (primitive since 
mgcd(M) = 1). Over Z3, the quadratic form reduces to 2q;^, which can represent and 2, 
but never 1, so that M and C are not conjugate over SL(2,Z3). This example highlights that 
the approach of [5j is of limited value for our problem. 

The following result is obvious from Proposition HI because Mat(2, Z„)^ -conjugacy induces 
a graph isomorphism in the sense of Definition [H compare Fact [71 

Corollary 2. Let M, M' be two integer matrices with the same trace and determinant, whence 
they share the same companion matrix. If, in addition, mgcd(M') = mgcd(M) = 1, the toral 
endomorphisms defined by M and M' have the same local statistics on all lattices L„. □ 

Example 6. Let us briefiy return to the matrices M, M^^ and C of Example [H They all 
have mgcd = 1, but are in different GL(2,Z) conjugacy classes. Nevertheless, Corollary [2l 
shows they all give the same local statistics, on all lattices L„. The same conclusion applies 
to the matrices M and M' of Example [3l 

Let us proceed to the general case when mgcd(M) = r > 1 from the point of view of 
local conjugacies. Consider the matrix M from (j20p . and let r = mgcd(M). We may now 
decompose M as 



(29) M = al + rN with N 



b 

c {d — a) 
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When r = 0, we have the trivial case M = al, which we now put aside by assuming r G N, 
in hue with Fact El so that mgcd(A^) = 1 by construction. 

Proposition 5. Let M = ("^) he an integer matrix with mgcd(M) = r 7^ 0. Then, for 
all integers n>2, and when taking all matrices mod n, M is Mat(2, Zn)^ -conjugate to the 
integer matrix ^ " 

Proof. By construction, the matrix N in the decomposition (j29p is still an integer matrix, 
with mgcd(A^) = gcd(6, c, (ti — a) ) = 1. By Proposition [H for all integers n > 2, is then 
Mat (2, Z„)^ -conjugate to its companion matrix 



Cat 



where both N and are taken mod n. 



be 

1 {d-a 



The claim now follows from the observation that this conjugacy extends to the one claimed 
by the structure of the decomposition M = al + rN. Indeed, the relation = ACj^A~^ 
together with (j29]) immediately implies AMA^^ = al + rC^. □ 

It is worth recalling that Lemma 1] does not apply to 2x 2-matrices, meaning that we 
do not have global conjugacy in general. The benefit of Proposition [5] is that we obtain a local 
conjugacy instead, and to a matrix with a well-defined element in the lower left corner. For 
the further arguments, we need a technical result, formulated within Mat(2,Z), hence prior 
to looking at the reductions mod n. 

Lemma 3. Let M,M' E Mat(2,Z) be the two matrices of (j20p and assume that det(M) = 
det(M') = D and tr(M) = tr(M') = T. Assume, in addition, that mgcd(M) = mgcd(M') = 
r G N. Then, r divides d — d' . 

Proof. When d = d' , the statement is trivially true, so we may assume that d ^ d' . Since 
r = gcd{b,c,d — a) = gcd{b' , c' , d' — a'), one has r^|6c and r'^\b'c'. Observe that 

{a-d'){d-d') = ad+{d' -T)d' = ad-a'd' = D + bc-D-b'c' = bc-b'c', 

which implies that r'^\[a — d') {d — d'). With m := d — d' ^ and d — a = kr, where k £ Zhy 
assumption, we now have to consider 

(30) I m{m — kr) 

and to show that this implies r\m. When m = kr, this is again clear, so assume m — kr ^ 0. 
Then, let = p^^^ ■ . . . ■ j?^'^ be the unique prime decomposition of into powers of distinct 
primes. Let ti be the highest power so that p*' \ m, which is a non-negative integer. 

Assume that ti < Si for some index i. By (f30]l . this implies \{m—kr) and thus, as pp\r 

and 2s j — ti > Si, also p^^m, in contradiction to the assumption that ti < Si. Consequently, 
we must indeed have ti > Si for all 1 < i < £. Since r = niP?' ^^^^ means that r|m as 
claimed. □ 



Proposition 6. Let M, M' be the matrices from (j20|) and assume that they have the same 
trace and determinant. Assume further that mgcd(M) = mgcd(M') = r. Then, for an 
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arbitrary integer n > 2, the reductions mod n of the matrices M and M' are Mat(2,Z„)^- 
conjugate. 

Proof. When r = 0, Fact [8] implies M = M' . So, let us assume r G N. By Proposition [H 
we know that matrices M and M' are Mat(2, Z„)^-conjugate to integer matrices N and A^', 
where 



in the sense that their reductions mod n satisfy the corresponding conjugacies. Our claim 
follows if we can show that also N and N' are Mat(2, Z„)^-conjugate in this sense. 



Consider the unimodular matrix A = ~r~ j , which is an integer matrix by Lemma [3] 
and hence an element of GL(2,Z). Using tr(M) = tr(M') and det(M) = det(M'), it is easy 
to check that AN = N' A holds, hence N' = ANA''^ within GL(2,Z). By Fact El this implies 



Let us investigate the local conjugacies a bit further, aiming at a converse of Proposition [H 

Proposition 7. Let M,M' E Mat(2,Z) be two integer matrices whose reductions mod n are 
Mat{2,Zn)^ -conjugate for some n > 2. Then, M and M' have the same determinants and 
traces mod n, and their matrix gcds r,r' generate the same ideal in Z„, meaning rZn = r"Zn- 

Proof. The statement about the determinants and traces is clear. For the claim about the 
matrix gcd, we can again use the idea of the proof of Lemma O up to the point where we 
conclude that r\r' and r'\r, now seen as divisibility properties within Z„. But k\i means that 
the principal ideal (i) = iZn is contained in (/c), so that we obtain (r) C (r') and (r') C (r), 
hence (r) = (r'). □ 

At this stage, one can formulate the following central result. 

Theorem 2. For two integer matrices M,M' £ Mat(2,Z), the following statements are 
equivalent: 

(a) The reductions mod n of M and M' are Ma.t{2, Zn)^ -conjugate for all n > 2; 

(b) M and M' satisfy the three relations det(M) = det(M'), tr(M) = tr(M') and 
mgcd(M) = mgcd(M'). 

Proof. The direction (b) =^ (a) follows directly from Proposition [6l 

For the converse direction, we may conclude from Proposition [7] that det(M) = det(M') 
mod n and tr(M) = tr(M') mod n for all n > 2. Consequently, we must have det(M) = 
det(M') and tr(M) = tr(M') (recall that k = i mod n G N means that A; — ^ is divisible by 
n, which simultaneously holds for all n G N only when k — £ = 0). 

For the third identity, let r = mgcd(M) and r' = mgcd(M') and assume that rZ„ = r'Z„ 
for all n > 2, but r ^ r' . Consequently, there is a prime p with r = q and r' = p* q\ 
t ^ s, such that g and g' are both relatively prime with p. Without loss of generality, we 
may assume that t > s, and then choose n = p^ . Clearly, both g and g' are then units 
in Z„ by construction. With (r) := rZ„, we then obtain (r) = {p^ g) = (p*) = (0), while 
(r') = {p^g') = (p*) / (0), in contradiction to the original assumption. □ 




and 





the local conjugacy claimed. 



□ 



20 



MICHAEL BAAKE, JOHN A. G. ROBERTS, AND ALFRED WEISS 



This theorem permits the following answer to the question for the local statistics of an 
integer matrix. 

Corollary 3. The complete local statistics of an integer matrix M , in the sense of Defi- 
nitionUl only depends on the three invariants det(M), tr(M) and mgcd(M). Two integer 
matrices with the same triple of invariants thus have the same local statistics on all lattices 
Ln- In particular, they have the same fixed point counts, both locally and globally. □ 

The result of Theorem [2] permits an interpretation in terms of Z, which is the inverse (or 
projective) limit of the rings Z„ over the positive integers ordered by divisibility, written as 
Z = lim Zn- It is also known as the Priifer ring, see [30] for details. 

Corollary 4. In the situation of Theorem [3 any of the two conditions is equivalent to the 
matrices M and M' being GL (2, Z)- conjugate. 

Proof. Clearly, a GL(2, Z)-conjugacy implies condition (a) of Theorem [2j Conversely, assum- 
ing (a), consider the inverse system of Mat(2, Z„) ^-subsets X(n) defined by 

X{n) = {A£ Mat(2,Z„)'' | AM = M'A} 

with n S N and ordered inductively by divisibility. All X{n) are non-empty by assumption 
and finite. Let X = limX(n) be the inverse limit, which is then non-empty as well. Any 
element in X achieves the required conjugacy via the corresponding projections. □ 

We conclude by noting: 

Remark 11. One consequence of Theorem [2] for SL(2,Z) matrices, noting p3p . is that such 
a matrix is conjugate to its inverse on all lattices individually, despite the fact that this is 
generally untrue on the torus, as in Example [H This relates to the investigations of quantum 
cat maps and their perturbations in [25]. 

Remark 12. We have described the three invariants for complete Mat(2, Z„) ^-conjugacy, 
whereas [3j presented the analogous result for GL(2, Z)-conjugacy. The conjugacy of M and 
M' in the special linear group over the p-adic integers is a related question, addressed in [5], 
though Example [5] above indicates the difference to GL(2, Z„)-conjugacy. 
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